A0-A178  417  MODIFICATION  OF  SEPARATING  FLOU(U)  OLD  DOMINION  UNTV 
NORFOLK  VA  DEPT  OF  OCEANOGRAPHV  C  E  GROSCH  FEB  87 
N88814-84-K-8585 


UNCLASSIFIED 


F/G  28/4 


DEPARTMENT  OF  OCEANOGRAPHY 
COLLEGE  OF  SCIENCES 
OLD  DOMINION  UNIVERSITY 
NORFOLK,  VIRGINIA  23508 


MODIFICATION  OF  SEPARATING  FLOW 


$ 


dtfc 

electei 

Mar  1  8  198/J 


Ci/7  D 


D 


Chester  E.  Grosch,  Principal  Investigator 


Interim  Summary  Report 

For  the  period  ended  June  30,  1985 


Prepared  for  the 
Department  of  Navy 
Office  of  Naval  Research 
Code  N00014 

800  North  Quincy  Street 
Arlington,  Virginia  22217-5000 


Under 

Research  Contract  N00014-84-K-0505 
Dr.  Robert  E.  Whitehead 

Scientific  Program  Officer  _ _ 

~Thu  dccurnont  he-:  bt-er>  o:«r.rcv«d 
loi  public  rtleost-  r-. H  >lo. 
distribution  is  t.nlbm-  H. 

Submitted  by  the  1 - - - - 

Old  Dominion  University  Research  Foundation 
P.  0.  Box  6369 
Norfolk,  Virginia  23508 


February  1987 


»7  3 


l  :*■  » •'  1 


MODIFICATION  OF  SEPARATING  FLOW 
By 

Chester  E.  Grosch* 

\  SUMMARY 

V 

The  usual  physical  description  of  separation  is  that  an  adverse~pres- 
sure  gradient  causes  the  deceleration  of  the  low  momentum  flow  near  the 
surface.  The  result  is  a  reduction  of  the  velocity  near  the  surface  and  a 
corresponding  decrease  in  the  shear  stress  at  the  wall  as  the  slope  of  the 
velocity  profile  decreases.  Finally,  at  some  point,  the  shear  stress  at  the 
wall  vanishes  and  the  flow  separates.  Downstream  of  this  point  the  flow 
near  the  boundary  reverses  and  a  "dividing41  streamline  breaks  away  from  the 
surface.  Below  this  streamline  the  net  flow  is  towards  the  separation  point 
and  above  this  streamline  it  is  away  from  the  separation  point.  This 
general  type  of  flow  field  can  generally  be  observed  near  the  rear  of 
bodies,  in  particular,  flow  separation  occurs  on  the  rearward  part  of  bodies 
with  a  blunt  shape. 

A  number  of  numerical  calculations  have  been  carried  out  in  order  to 
study  the  effect  of  transverse  grooves  on  flow  separation  and  incipient  flow 
separation.  These  calculations  are  intended  to  throw  light  on  the  dynamics 
of  the  flow  toward  the  rear  of  two-  and  three-dimensional  bodies.  In  the 
most  general  case,  body  curvature  should  be  included,  and  it  would  be  neces¬ 
sary  to  calculate  the  flow  past  the  body  in  a  curvilinear  coordinate  system. 
This  would  be  a  massive  undertaking  for  general  geometries.  In  order  to 
isolate  the  effect  of  the  grooves  on  flow  separation,  apart  from  the  effects 
of  body  curvature,  it  was  decided  to  use  planar  geometry  only. 


♦Professor,  Department  of  Oceanography,  Old  Dominion  University,  Norfolk, 
Virginia  23508 


That  is,  the  problem  was  taken  to  be  the  flow  of  an  incompressible  fluid  in 
a  two  dimensional  geometry  past  a  plane  wall.  There  is  an  externally  im¬ 
posed  pressure  gradient  and  a  slot  in  the  boundary.  Calculations  were  car¬ 
ried  out  with  the  slot  closed  as  well  as  open.  The  equations  solved  are  the 
two-dimensional  Navier- Stokes  equations  for  the  flow  of  an  incompressible 

fluid.  These  are  written  in  terms  of  the  velocity,  u  =  (u,  v),  and  vor- 
ticity,  c .  These  are 


u  +  v  =0, 
x  y  * 


v  -  u  =  c, 
x  y 


«y=Re->  (cxx*tyy). 


together  with  appropriate  boundary  conditions.  This  set  of  differential 
equations  is  replaced  by  a  corresponding  set  of  difference  equations.  The 
difference  equations  used  are  of  the  compact  form,  that  is  only  variables  in 
a  cell  or  on  the  boundaries  of  a  cell  are  needed.  This  formulation  has  some 
advantages  over  more  conventional  difference  schemes.  In  particular,  it  is 
trivial  to  use  nonuniform  grids  as  was  done  in  these  calculations.  The 
nonuniform  grids  were,  of  course,  used  to  give  better  resolution,  in  the 
near  wall  region  and  in  the  vicinity  of  the  slot,  than  elsewhere. 

Defining  the  difference,  5,  and  averaging,  u,  operators  by 


i,j  S  (fi  +  1/2,  j  _  fi  -  1/2,  j) 


;\/AX 


A  detailed  derivation  of  these  equations  is  given  elsewhere  (Gatski,  Grosch, 
and  Rose,  1982). 

The  geometry  of  the  flow  for  which  calculations  were  carried  out  was  a 
solid  boundary,  both  with  and  without  a  slot.  There  was  an  inflow  boundary 
and  an  outflow  boundary.  At  the  "top"  there  was  a  freestream  boundary. 

At  the  inflow  boundary  the  normal  component  of  the  velocity,  u,  and 

the  vorticity,  c,  are  given.  In  all  of  the  calculations  they  were  taken 
to  be  those  of  the  zero  pressure  gradient  Blasius  boundary  layer.  Calcula¬ 
tions  were  carried  out  for  only  one  inflow  Reynolds  number.  The  Reynolds 
number,  based  on  the  free  stream  speed  and  the  displacement  thickness  of  the 
boundary  layer  was  approximately  380.  This  was  chosen  so  that  the  basic 
flow  would  be  that  of  a  stable  laminar  boundary  layer.  On  the  solid  wall, 
and  on  the  walls  of  the  slot  when  it  was  taken  to  be  open,  the  boundary 
conditions  were  taken  to  be  that  the  normal  and  tangential  components  of  the 
fluid  velocity  was  zero.  On  the  outflow  boundary  approximate  boundary 
conditions 

dv/dt  =  dc/dt  =  0  ,  (16) 

were  used.  Finally  on  the  free  stream  boundary  on  the  "top"  the  tangential 
component,  u,  of  tne  velocity  was  specified  along  with  an  outflow 
condition 

dc/dt  =  0  ,  (17) 

on  the  vorticity.  Specifying  u  as  a  function  of  the  downstream  distance, 
x,  is  equivalent  to  specifying  the  pressure  gradient  3p/3x.  This  is  clear 
if  one  notes  that  in  the  free  stream 


(18) 


p+ip(u2+v2)=  constant. 
2 


thus  the  specification  of  u(x),  and  thus  v(x)  through  v.u  =  0, 

k  •  Vx  u  =  c,  effectively  specifies  p(x).  In  fact,  noting  that 
|v|/|u|  «1, 

-  3  p/3x  -  I  3(u2)/3x.  (19) 

P  2 


With  these  boundary  conditions  the  flow  is  completely  specified. 
However,  there  was  one  difficulty.  If  u(x)  on  the  "top"  was  given  so  that 
3 p/3x  was  unfavorable  and  separation  occurred,  one  could  not  specify 
"outflow"  boundary  conditions  on  the  downstream  boundary  because  it  was,  in 
part,  an  inflow  boundary.  The  proper  inflow  conditions  on  these  sort  of 
boundaries  are  not  known,  even  in  principle.  In  order  to  circunvent  this 
problem,  u(x)  or  the  "top"  boundary  was  chosen  so  that  the  flow  had  a 
region  of  zero  pressure  gradient,  followed  by  an  unfavorable  pressure 
gradient,  then  a  region  of  zero  pressure  gradient.  This  resulted  in  the 
formation  of  a  region  of  decelerating  and  near  separating  flow  or  the 
formation  of  a  separation  bubble  on  the  boundary.  In  all  cases,  the  flow  at 
the  outflow-boundary  was  that  of  an  attached  flow. 

The  solution  algorithm  for  the  set  of  difference  equations  is  described 
in  some  detail  by  Gatski,  et  al .  (1982)  and  will  not  be  discussed  further 
here.  It  will  suffice  to  say  that  with  the  equations  and  boundary 
conditions  described  above,  steady  state  solutions  were  obtained  in  all 
cases.  This  was  true  whether  or  not  the  boundary  slot  was  open  or  closed. 


However,  in  some  cases  many  time  steps  were  required  in  order  to  reach  a 
steady  state. 

Calculations  were  carried  out  for  a  number  of  different  cases.  In  all 

of  these  u(x)  =  1.0  at  inflow  and  for  some  distance  thereafter.  It  is 

then  reduced  to  a  u^-p  and  held  at  that  value  up  to  the  outflow. 

Calculations  have  been  carried  out  for  a  range  of  values  of  u  .  =  1.0, 

mm 

to  0.70.  This  gave  a  range  of  adverse  pressure  gradients.  In  all  cases  the 
"opening"  of  a  slot  in  the  boundary  significantly  modified  the  flow  in  the 
neighborhood  of  the  slot.  If  the  slot  were  sufficiently  deep,  say  greater 
than  one  boundary  layer  thickness,  there  was  no  significant  increase  in  the 
drag  due  to  the  slot.  However,  the  flow  was  locally,  accelerated  due  to  the 
presence  of  the  slot  and  this  had  an  important  effect  on  tne  separating  or 
near  separating  flow. 
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